Search for: All records

We provide an approach to sample rare events during classical ab initio molecular dynamics and quantum wavepacket dynamics. For classical AIMD, a set of fictitious degrees of freedom are introduced that may harmonically interact with the electronic and nuclear degrees of freedom to steer the dynamics in a conservative fashion toward energetically forbidden regions. A similar approach when introduced for quantum wavepacket dynamics has the effect of biasing the trajectory of the wavepacket centroid toward the regions of the potential surface that are difficult to sample. The approach is demonstrated for a phenol-amine system, which is a prototypical problem for condensed phase-proton transfer, and for model potentials undergoing wavepacket dynamics. In all cases, the approach yields trajectories that conserve energy while sampling rare events. 

Abstract The time‐dependent Schrödinger equation can be rewritten so that its interpretation is no longer probabilistic. Two well‐known and related reformulations are Bohmian mechanics and quantum hydrodynamics. In these formulations, quantum particles follow real, deterministic trajectories influenced by a quantum force. Generally, trajectory methods are not applied to electronic structure calculations as they predict that the electrons in a ground‐state, real, molecular wavefunction are motionless. However, a spin‐dependent momentum can be recovered from the nonrelativistic limit of the Dirac equation. Therefore, we developed new, spin‐dependent equations of motion for the quantum hydrodynamics of electrons in molecular orbitals. The equations are based on a Lagrange multiplier, which constrains each electron to an isosurface of its molecular orbital, as required by the spin‐dependent momentum. Both the momentum and the Lagrange multiplier provide a unique perspective on the properties of electrons in molecules.